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^ ' Abstract. We derive the galaxy luminosity functions in V-, R-, I-, and K-bands of 

' the cluster EIS 0048 at z ~ 0.64 from data taken at the ESO Very Large Telescope. 

' The data span the restframe wavelength range from UV, which is sensitive to even 
low rates of star formation, to the NIR, which maps the bulk of the stellar mass. 
, By comparing our data and previous results with pure luminosity evolution 

(-H , models, we conclude that bright {M < M* + 1) cluster galaxies are already as- 
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sembled at z ~ 1 and that star formation is almost completed at a ~ 1.5. 
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1- Introduction 

One of the main issues in theoretical and observational research is to understand the 
physical processes driving the formation and evolution of bright massive galaxies in clus- 
ters, and to constrain the relative time scales. 

The evolution of the cluster galaxy sequence (the colour-magnitude relation) has been 
studied for large cluster samples up to redshift z ~ 1 (Aragon-Salamanca et al. 119931 
Lubin lTMSI Ellis et al. lTWTI Stanford et aL lTMSl hereafter SED98; Nelson et al. lMTTl 
hereafter NGZOl). The results are consistent with a monolithic collapse scenario (see e.g. 
Larson [l974|l in which galaxies form at high redshifts and subsequently undergo passive 
evolution. However, the evolution of colours only inform on the epoch when the bulk of 
stellar mass formed, while cluster bright galaxies [M < M* + 1, mostly early-type) could 
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have been also assembled recently (z < 1) from mergers of smaller units, at least as long 
as the merging processes do not induce strong star formation. 

A different approach consists in the study of the evolution of the cluster luminosity 
function (CLF). In particular, the near-infrared (NIR) CLF can be used to assess the 
assembly history of galaxies, because the NIR light mainly informs on galaxy mass. From 
the analysis of the values of M* in the K-band as a function of rcdshift for a sample of 
38 clusters in the range 0.1 < z < 1, De Propris et al. (jl999t hereafter DPSE99) found 
a trend consistent with passive luminosity evolution (PLE) (see also Nakata et al. 120011 
hereafter NKYOl). DPSE99 conclude that the mass function of bright cluster galaxies is 
invariant at z < 1, and that the assembly of those galaxies is largely complete by z ~ 1. 
Spectroscopic data also show that massive galaxies already exist at least up to z = 0.83 
(van Dokkum et al. lHMjl . 

In this work we derive the CLFs for EIS 0048 at z 0.64 in the V-, R-, I-, and K-band 
(UV to NIR restframe) . The cluster membership has been assessed with the photometric 
redshift technique up to M — M* ~ 1.5 — 2.5 (according to the different depth of each 
band). Since no other selection has been applied, the samples are not biased toward a 
particular galaxy population. 

The paper is organized as follows. In Sect. 2 we introduce the photometric data, 
discuss the background subtraction, the selection of cluster members, the completeness 
of the samples, and obtain the CLFs. In Sect. 3 we discuss the results in terms of galaxy 
formation and evolution. In Sect. 4 we give a summary of the paper. In this work we 
assume Hq = 70 km s^^ Mpc^^, = 0.3 and Ha = 0.7. 

2. Derivation of the luminosity functions 

The photometric observations of the cluster of galaxies EIS 0048 were carried out at the 
ESO Very Large Telescope (VLT) during two observing runs on August 2001. All the 
nights were photometric with excellent seeing conditions. 

The data include VRIK imaging taken with the F0RS2 and ISAAC instruments, 
respectively. The VRI images consist of a single pointing of 6.8x6.8 arcmin^ for each band, 
while for the K-band a mosaic of four pointings covers a total area of 4.9x4.9 arcmin^. 
Further details on the observations and on the data reduction can be found in La Barbera 
et al. plHHl hereafter LBMI03). 

2.1. Magnitudes and completeness 

Total magnitudes were computed by means of the software SExtractor (Bcrtin & Arnouts 
1996). For each object we obtained Kron magnitudes (mx) within an aperture of diameter 
a ■ rK, where ik is the Kron radius (Kron 1980). We chose a = 2.2, for which mK is 
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Fig. 1. Completeness limits of the photometric catalogues. Total magnitudes (magfot) are 
plotted against the magnitudes in the detection aperture (magce/;)- The vertical and hor- 
izontal lines correspond to the detection threshold and to the completeness magnitudes 
respectively. 

expected to enclose 92% of the total flux, and we computed the total magnitudes by 
adding 0.08 mag to mK (see LBMI03 for details). 

The formal completeness of the catalogues was estimated following the method of 
Garilli et al. 1)19991 see also Busarello et al. 2002). This method consists in the determi- 
nation of the magnitude at which the objects start to be lost since they are below the 
brightness threshold in the detection cell ^. To estimate the completeness limit, mag- 
nitudes in the detection cell were compared to the total magnitudes. The comparisons 
are shown in Fig. ^ where the vertical lines correspond to the detection threshold and 

^ In the present case, the detection cell is an aperture with diameter of 5 pixels. 
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Table 1. The completeness magnitudes, m'^ is the completeness magnitude as estimated 
on the photometric catalogues according to Garilli et al. H1999|l . m^p defines the central 
position of the faintest magnitude bin with completeness higher than 50% in the sample 
with photometric redshifts (see Sect 12. 2|) . The bin width is 0.5 mag. 
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the horizontal lines mark the completeness limits (taking into account the scatter of the 
relation). The completeness magnitudes are reported in Tabled 

2.2. Background subtraction 

The LFs of cluster galaxies at low and intermediate redshifts are usually obtained by 
statistical subtraction of the background and foreground contribution (as derived from 
a control field) to the galaxy counts in the cluster field. In this way all the photometric 
data up to the magnitude limit can be used. On the other hand, some problems become 
severe as the cluster redshift increases: a) cluster galaxies become fainter; b) the number 
of galaxies in the control field becomes higher; c) clusters are younger and therefore 
the number of cluster galaxies is smaller. The consequences are that the contrast of 
cluster galaxies over the field decreases as a function of redshift, and that the Poissonian 
fluctuations of the field as well as field-to-field variations become higher than the signal 
to be detected (see also NKYOl). 

To analyse this problem, we compare a known CLF translated to different redshifts, 
with the uncertainties on field galaxy counts. We use the LF of the rich cluster AC 118 
at z = 0.31, iVcLu("T'), in the K-band from Andreon H2(JU1|I . and the galaxy counts 
A'hdfs(™) and iVcDF("^) from two different control fields: the Hubble Deep Field South 
(HDFS) and the Chandra Deep Field (CDF) respectively, as obtained by Saracco et al. 

We define the cluster signal-to-noise ratio CSNR as: 

CSNRim) = ^^^^("^^ , (1) 

v/(A^HDFs(m) - A^cdf(to))2 + A^cLu(m) 

where A'clu("t.) is the signal we aim to measure and 



\/(A^HDFs(™) — A''cdf(w))^ + iVcLu(w) is an estimate of the uncertainties on the 
field galaxy counts plus the Poissonian errors on the CLF. 
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Fig. 2. The cluster signal-to-noise ratio (see text and Eq.QJ in a rich cluster at z = 0.31 
(filled dots and continuous line) and at z = 0.64 (open squares and dotted line). 



We computed the CSNR for AC 118 at z = 0.31 (filled circles and continuous line in 
Fig. [5J) and then shifted it to z = 0.64 (open squares and dotted line) according to the 
PLE model described in Sect. El As can be seen in Fig. |3 for z = 0.31 CSNR{m) > 2 
at aU magnitudes and CSNR{m) > 3 at to > M* + 0.5 {M* = 15.3, see Andreon lM7T|l . 
At z = 0.64, CSNRim) < 2 at all magnitudes, implying that the signal to be detected 
would be almost equal to the noise. 

An additional problem arises from two facts: a) AC 118 is a cluster richer than 
EIS0048; b) by shifting the LF of AC 118 at a higher redshift, we overestimated the ac- 
tual number of detected galaxies. As consequence, our estimate of CSNR{m) at z — 0.64 
must be considered as an upper limit, and adopting the procedure of statistically sub- 
tracting the background and foreground contribution from a control field would prevent 
an accurate estimate of the CLFs. It should also be noted that we do not have access to 
a control field in the immediate proximity of EIS 0048 in order to minimize uncertainties 
from field-to-field variations (as in NGZOl), and that we cannot derive a cumulative LF 
by combining the signals from different clusters at similar redshifts (as in DPSE99). 

Therefore, in order to increase the contrast of cluster galaxies over field galaxies, 
the following approach was adopted. We first selected the cluster members through the 
photometric redshift technique, as detailed in LBMI03. In this way it is possible to isolate 
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Fig. 3. The completeness as a function of the magnitude of the samples used to derive 
the photometric redshift distribution of the galaxies in the EIS 0048 field. In each panel 
the vertical dotted line indicates the magnitude corresponding to 50% completeness. The 
CLFs are computed up to the faintest magnitude bin brighter than these limits. 



galaxies in a narrow interval AziO.l around the cluster redshift, thus greatly decreasing 
the noise produced by field galaxies, while leaving the signal untouched. 

To achieve a reasonable accuracy in the redshift estimate, we considered only galaxies 
with signal-to-noise ratio SjN > 5 in at least three bands. For each band, we computed 
the completeness function of the sample as the ratio of the number of galaxies with 
SjN > 5 in at least three bands to the total number of galaxies in the catalog, as a 
function of the magnitude (see Fig. |3J|. The CLFs have been computed in each band by 
considering all the galaxies with photometric redshift up to the last magnitude bin with 
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Fig. 4. The CLFs of EIS 0048 from V- to K-band (filled dots). In each panel the best-fit 
Schechter function is also shown (solid line). 

completeness higher than 50%, as reported in Tabled Galaxy counts have been corrected 
for the estimated incompleteness. 

After the previous selection, field contamination is still present in the peak of the 
redshift distribution around z ~ 0.64. To account for the remaining contamination, we 
took advantage of the VIRMOS preparatory photometric survey as control sample (see 
LBMI03). By selecting galaxies in the redshift range z £ [0.54,0.74] (LBMI03) in the 
VIRMOS control field, we obtained the luminosity distribution of contaminant galaxies, 
and subtracted it from galaxy counts in the same redshift range in the cluster field. At 
K < 20 the ratio of cluster to field galaxies without any photometric redshift selection 
is ~ 0.5 (from Saracco et al. 120011 data in the CDF), whereas the same ratio after 
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Table 2. The values of M* computed from the fits of the Schechter function to the CLFs 
of EIS 0048. M* is obtained by fixing a = -0.9. 



Band M* 
mag 

V 22.70 ±0.25 
R 21.27 ±0.20 
I 20.24 ±0.18 
K 17.18 ±0.23 



the photometric redshift selection is ^ 6.6 (from the VIRMOS data), thus proving the 
effectiveness of our procedure. 

2.3. Luminosity functions 

We modelled the CLFs with a weighted parametric fit of the Schechter function, fixing 
the faint end slope to have a — —0.9. We chose not to fit a because the samples with 

photometric redshifts are only complete to M* + 1.5 2.b (according to the band). 

Moreover, a = —0.9 is the value measured in the NIR for the Coma cluster &i K < K* + ?> 
(De Propris et al. ll998jl and is also the value adopted by DPSE99 to study the behaviour 
of K* as a function of the redshift. The errors on the CLFs were computed by taking into 
account Poissonian fluctuations of cluster galaxies counts and of the background counts 
in the redshift range z e [0.54,0.74]. In Fig. 0] we show the CLFs from V- to K-band 
(filled circles) and the best-fitting Schechter functions (solid lines). The values of M* 
computed by the fits are listed in Table HJ 

In the K-band we obtain M* — 17.18 ± 0.23, in agreement within the errors with the 
result of DPSE99 at the same redshift {K* = 17.57 ± 0.41). By comparing the values of 
M* in I- and K-band, we obtain /* — K* = 3.06 ± 0.28 for the bright cluster galaxies at 
z — 0.64, in agreement with the findings of SED98 and NGZOl at similar redshifts. 

3. Luminosity evolution of bright (M < M* + 1) cluster galaxies 

We discuss the luminosity evolution of bright galaxies in clusters on the basis of 
the behaviour of the V- and K-band characteristic magnitudes M* with redshift. 
Complementary data for the K-band arc taken from DPSE99, Andreon l2001l and NKYOl, 
while for the V-band are taken from Busarello et al. 120021 The different data sets are 
represented in Fig.^with different symbols. The points at z = 0.31 in the K- and V-band 
are obtained from the luminosity functions of Andreon H20()l|l and Busarello et al. I)2()02|l . 
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Fig. 5. Evolution of the Schechter parameter M* in K- and V-bands as a function of 
redshift. In the K-band the data come from DPSE99 (empty circles), Andreon 1)20011 
empty square), NKYOl (empty pentagon), and from this work (filled circle). In the V- 
band the data come from Busarello et al. H2002L empty triangle) and from this work 
(filled circle). In the upper panels, the dashed lines represent the NE model. The PLE 
models (solid lines) at different formation redshift Zf are shown in decreasing order of 
Zf from top to bottom. In the lower panels the PLE models (dotted lines) are shown 
as a function of the redshift zgg for which where galaxies form the 99% of their stars 
(decreasing from top to bottom). 

assuming a — —0.9 up to M ^ M* + 2. In this way all the numerical estimates of M* 
have been obtained with the same value of a and can be straightly compared. 

In Fig. El no-evolution (NE, pure /c-correction) and PLE models are also shown. The 
galaxy models are constructed using Bruzual & Chariot H199:^|l templates with a Scalo 
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initial mass function, solar metallicity, and an exponential star formation law e^*/"^. We 
assumed the time scale t = 1 Gyr, which is suitable to reproduce the colour evolution of 
early type galaxies (see Bruzual & Charlot ri993l Merluzzi et al l2U03fl . The NE predictions 
use model spectra with formation redshift Zf — 5 (the redshift of the onset of SF). 
The predictions from PLE are given for different Zf. All the models are normalized to 
K* = 10.90 and V* = 14.17 at the redshift of Coma (Bower et al. ITM^ De Propris et 

ai.CMHI). 

As widely discussed in literature (e.g. DPSE99; NKYOl), the behaviour of K* is 
inconsistent with a NE model. Our data points at z = 0.64 confirm this result (at 
~ 2.50-), and extend it in the V-band (at ^ Aa). The evolution in the K-band can be 
described by PLE models at least up to z 1 (the point at z = 1.2 in Fig. [S] should 
be taken with caution since NKYOl were not able to estimate the field contamination). 
Since the K-band samples the bulk of stellar mass, this result suggests that the assembly 
of bright galaxies in clusters is complete at z ^ 1 (DPSE99). 

The K-band data alone cannot discriminate among models with different formation 
redshifts. At z = 0.64 the V-band samples the UV restframe and is therefore sensitive 
to even low levels of star formation activity. As can be seen in the upper right panel in 
Fig- El the evolution of V* is in agreement with PLE models with zy > 3. This result 
depends on the time scale of star-formation history, but we can use our point at z — 0.64 
in the V-band to constrain the time when galaxies ceased to form stars. This is possible 
because the magnitude evolution subsequent to this time is almost independent from the 
preceding star-formation history. 

In the lower panels of Fig. we plot the PLE models corresponding to the different 
epochs at which galaxies formed 99% of their present stars. As it is apparent from the 
Figure, the V-band data constrain the star formation in cluster bright galaxies to end at 
z > 1.5. 

4. Summary and conclusions 

We derived the galaxy LF in V-, R-, I-, and K-bands of the cluster EIS 0048 at z ^ 0.64 
from new photometric observations carried out at ESO VLT with F0RS2 and ISAAC. 
Cluster members have been selected with the photometric redshift technique in order to 
enhance the contrast among cluster- and background-foreground galaxies. The remaining 
field contamination has been estimated as detailed in LBMI03. 

The CLFs have been obtained up to M — M* ^ 1.5 — 2.5 (according to the waveband). 
We modelled the CLFs with a weighted parametric fit of the Schechter function, fixing 
the faint end slope to be a = —0.9. The values of M* obtained in I- and K-bands are 
in agreement with the analysis of DPSE99 and NGZOl at similar redshifts. We collected 
results from literature and introduced PLE models with different formation redshifts in 
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order to discuss the evolution of M* in the V- and K-bands. The key factors driving 
the evolution of bright galaxies in the scenario of hierarchical structure formation are the 
epoch when the bulk of stellar populations is formed, the cosmological time when mergers 
are effective to assemble the galaxies, the amount of star formation induced by mergers, 
and the age of the youngest stars. Since the evolution of K* is consistent with a PLE 
scenario at least up to z ^ 1, we can conclude that mergers must have already assembled 
bright (M < M* + 1) cluster galaxies at this redshift (see also DPSE99 for a thorough 
discussion). At z = 0.64 the V-band samples the UV restframe wavelength region and 
is sensitive to even low levels of SF. We find that the SF in bright cluster galaxies has 
to be almost completed at z 1.5, whereas the formation redshift is zy > 3 assuming 
r ~ 1 Gyr as the time scale for SF. These results lead to conclude that the structure 
and the stars of bright cluster galaxies must have been formed between z = 4 ± 1 and 
z = 1.2 ±0.2. 
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